Fully gapped superconductivity in SrNi 2 P2 
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We investigated the superconducting gap structure of SrNi2P2 (T c = 1.4 K) via low-temperature 
magneto-thermal conductivity re(T, H) measurements. Zero field thermal conductivity re decreases 
exponentially re oc exp(—aT c /T) with a = 1.5, in accord with the BCS theory, and rolls over to 
a phonon-like re oc T 3 behavior at low temperature, similar to a number of conventional s-wave 
superconductors. In addition, we observed a "concave" field dependence of the residual linear term 
ko(H)/T. These facts strongly rule out the presence of nodes in the superconducting energy gap 
of SrNi2P2- Together with a fully gapped Fermi surface in the superconducting state of BaNi2As2 
(T c = 0.6 - 0.7 K) , demonstrated in our recent work, these results lead us to stipulate that fully gapped 
superconductivity is likely to be a universal feature of Ni-based pnictide superconductors. 
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I. INTRODUCTION 

Since the discovery of superconductivity in 
LaFeAs(0,F)fi one of the so-called Fe-pnictides, a 
considerable number of experimental and theoreti- 
cal investigations have been performed on this class 
of materials, which includes compounds with high 
superconducting transition temperatures (T c ) up to 
^56Kf2r— in proximity to magnetism accompanied 
by structural transitions, and in a variety of crystal 
structures which share FePn layers (Pn = As, P, Sc 
etc.)£ The important issue of the superconducting 
pairing symmetry in Fe-pnictides is still controversial. 
Conclusions for the gap structure range from a single 
s-wave-like gap^ to multiple gaps^ to the existence 
of nodes in the gap, depending not only on the probe 
used to make the measurement, but also on a particular 
sample composition^— Consequently, the gap function 
in Fe-pnictides is proposed to be dependent on the 
structure of the Fermi surface and to be non-universal* 1 - 1 - 
Theoretically, a sign-reversing s± model, based on the 
unique Fermi surfaces in Fe-pnictides ji 2 ^— was proposed 
by several groups to resolve many apparent experimental 
discrepancies^ Investigations of the related families 
of compounds may help to resolve the issue of the 
superconducting gap symmetry and the origin of the 
pairing mechanism in Fe-pnictides. 

Ni-analogs (Ni-pnictides), with the same structure as 
Fe-pnictides, also superconductJ^r— In addition, some 
Ni-pnictides share several properties with Fe-pnictides, 
such as structural transitions from tetragonal to a lower 
symmetry crystal lattice ; 24 ' 25 enhancement of T c by 
doping ; 26 ' 27 and the importance of low-dimensionality 
for achieving higher T c i 25 ' 28 On the other hand, there 
are crucial differences including the magnitude of T c , 
which does not exceed 5K in any of the Ni-based com- 
pounds, the absence of magnetic ordering,— and a more 
three dimensional structure of the Fermi surface in Ni- 
pnictides j2£ Correlating the commonalities and, more im- 



portantly, the differences in physical properties between 
Fe- and Ni-pnictides could hold important clues to, in 
particular, the origin of the high T c in Fe-pnictides. 

In our recent thermal conductivity k and specific 
heat C studies of BaNi 2 As 2 ,22 with T C = 0.7K and 
a structural transition To=130K;2i we have estab- 
lished fully gapped superconductivity in this Ni-pnictidc. 
Similar conclusions were recently arrived at for La- 
NiAs(0,F) via nuclear quadrupolc resonance (NQR) 
measurements^ While the electronic structure for vari- 
ous Ni-pnictides are qualitatively similar ) 29 ' 32 " — the de- 
tailed differences between them allow us to check for 
universality of the superconducting properties. Here, 
we report a thermal conductivity study of another 
Ni-pnictide, SrNi 2 P 2 (T = 325K, T C = 1.4K)^ Low- 
temperature thermal transport is an established pow- 
erful tool for investigating superconducting properties. 
In particular, the low-temperature magnetic field de- 
pendence of the thermal conductivity k is instructive 
to detect whether the Fermi surface is fully gapped in 
the superconducting state 30 i 35 ~ iS or there are nodes in 
the gap i 41 ' 42 Thermal conductivity of SrNi 2 P 2 reported 
here is qualitatively similar to that of BaNi 2 As 2 , which 
leads us to suggest that fully gapped superconductivity 
might be a universal feature in Ni-pnictides, as opposed 
to a large variation in a gap structure observed in Fe- 
pnictides. 



II. EXPERIMENTAL DETAILS 

SrNi 2 P 2 forms in the ThCr 2 Si 2 -type tetragonal struc- 
ture. Single crystal samples were grown in Sn flux 
as described in previous reports! 25 ' 43 Thermal conduc- 
tivity measurements were performed via standard one- 
heater-two-thermometers technique on a plate-like crys- 
tal (~ 2 x 1 x 0.3 mm 3 ) in a range of 50 mK to 3K using 
an S.H.E. corp.'s dilution refrigerator equipped with a 
9T superconducting magnet. The heat current q was 
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FIG. 1: (Color online) Main figure: Temperature dependence 
of the thermal conductivity re(T) of SrNi2P2 in zero field for a 
heat current q || [100], on a double logarithmic plot. An arrow 
indicates T c (=1.4K), determined from the resistivity and 
specific heat data. A dotted curve shows a fit to k(T) with 
a BCS curve defined as k = Cexp(— aT c /T). A dashed-dotted 
line represents the calculated upper limit of the phonon con- 
ductivity Kpjj (ocT 3 ) based on the scattering off the sample 
boundary. Inset: measured (open symbols) and estimated 
electronic thermal conductivity in normal state (solid trian- 
gles). A straight line is a guide to the eye for kocT. Applied 
field H || [010], and the heat current q || [100]. 



applied along [100], the longest dimension of the sam- 
ple. Pt wires, spot-welded to the sample, provided ther- 
mal links to the heater, thermometers, and the bath. 
Meanwhile, superconducting NbTi wires provided elec- 
trical connection to the heater and Ru02 thermometers, 
while thermally isolating them from the support frame. 
Electrical resistivity was measured down to 0.4 K in a 
physical property measurement system (PPMS: Quan- 
tum Design), using the crystal (with the same electrical 
contacts) that was used in thermal conductivity measure- 
ments. Magnetic field was always applied along the sec- 
ond longest dimension of the sample (within the plate), 
H || [010], resulting in a H _L q arrangement. Measure- 
ments were performed upon field cooling from above T c 
to achieve a uniform magnetic flux distribution. Mag- 
netic field history dependence was explored in a separate 
set of experiments described below. 



III. RESULTS 

Figure [TJ shows the temperature dependence of ther- 
mal conductivity re(T) of SrNi2P2 in zero field for 



the heat current g||[100]. In the normal state, as 
shown in the inset of Fig. [TJ re(T) follows an approx- 
imately T-linear variation both above T c in zero field 
and in the whole temperature range measured in a 
field of 0.1 T (>iJ c2 (0)=0.039T^). Solid triangles 
represent the electronic conductivity re in the normal 
state, estimated from resistivity via the Wiedemann- 
Franz law: K e = LoT/p, with the Lorenz number 
L = 2.44 x 10~ 8 Wf2/K 2 . Good agreement between the 
measured re and the estimated n c in the normal state 
in fields above H C 2 implies that the heat transport in 
the normal state of SrNi2P2 is dominated by the elec- 
trons, and the phonon conductivity is negligible. In the 
superconducting state below T c , re(T) follows an expo- 
nential form (~ exp(— aT c /T), expected from the BCS 
theory, with a = 1.5 down to 0.5 K, as indicated by the 
dotted curve. Below 0.5 K, re(T) starts to deviate from 
the exponential T-dependence and, below 0.3 K, it fol- 
lows a power law T a dependence, with a = 3. A reocT a 
behavior, with 2 < a < 3, has been reported in a 
number of conventional s-wave superconductors, and is 
commonly attributed to a dominant phonon conductiv- 
ity: K ph- K p h can overcome the exponentially- reduced re e 
at low-temperature, even when re p h <C re e i n the normal 
state, as in SrNi2P2- The upper limit of the phonon 
thermal conductivity K^ h (oc T 3 ), is determined by the 
scattering off the sample boundaries.— In fact, using 
the coefficient /3 = 4.6J/K 4 m 3 of the low temperature 
T 3 term of the phonon specific heat, mean phonon ve- 
locity, (v) = 3 x 10 3 m/s— , and the phonon mean free 
path =\/4a&/7r = 6.3 x 10 -4 m, with the dimensions 
a = 1.2 mm and b = 0.3 mm of the sample cross-section, 
we obtain re^ h = ^T 3 (v)l ph = 2.9 x T 3 W/K 4 m, shown 
by a dashed-dotted line in the main panel of Fig. 1. This 
estimate is large enough to account for the T 3 variation 
of the low-temperature conductivity, and the discrepancy 
between the estimate and experimental data indicates 
an additional scattering mechanism, which reduces the 
phonon mean free path below the value above determined 
by the sample cross-section. 

The thermal conductivity of SrNi2?2 displays features 
common to other conventional s-wave superconductors 
and BaNi2As2. In figure [2] we compare re vs T/T c 
in zero field for SrNi 2 P 2 , BaNi 2 As 2 (#l, T C = 0.68K),^ 
BaNi 2 As 2 (#2, T C = 0.62K)^ and the s-wave supercon- 
ductor AL 44 ' 46 Despite widely varying absolute values, 
the overall temperature variation of re(T) in SrNi2P2 
and in another Ni-pnictidc BaNi 2 As 2 resembles that of 
a number of conventional s-wave superconductors, e.g. 
Al, in terms of (i) a T-linear variation due to electrons 
in the normal state, (ii) a BCS-like exp(— aT c /T) behav- 
ior with a ~ 1.3-1.5 in the superconducting state be- 
low T C ) 44 ' 47 and (iii) a T Q -dependence with cv = 2-3 in 
the lowest temperature region. In contrast, re(T) of Fe- 
pnictides exhibits a distinct rise when the sample enters 
the superconducting state4^r— The rise at T c in re(T) in 
Fe-pnictides superconductors is reminiscent of thermal 
conductivity in high-T c cuprates^ and in heavy fermion 
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FIG. 2: (Color online) Comparison of zero-field k vs T/T c 
in SrNi 2 P 2 , BaNi 2 As 2 (#1, T C =0.68K),^ BaNi 2 As 2 (#2, 
T c — 0.62 K),— and a conventional s-wave superconductor Al 
(T c = 1.2 K)4£ For Al, k is multiplied by 0.02. Solid lines are 
guides-to-the-eye representing ocT" (n= 1,2,3). 



compounds, such as CeCoIn 5 ^ which is a result of a re- 
markable increase of the electronic mean free path in the 
superconducting state that overwhelms the reduction in 
the density of states. 

Next, we turn to the magnetic field dependence of 
the low-temperature thermal conductivity in SrNi2P2- 
The magnetic field dependence of the residual linear 
term ^-/^\ T ^ has been used in the past (e.g. in 
the case of BaNi2As2) to identify the structure of the 
superconducting gap symmetry. Namely, a concave 
(oc ^/~HeyLj>(—b\jH C 2/H) dependence of the residual lin- 
ear term at low magnetic fields indicates a fully gapped 
Fermi surface in the superconducting state, while a con- 
vex (oc y/~H) field dependence indicates the presence of 
nodc(s) in the energy gap. Figure E3a) shows k/T 
vs T of SrNi 2 P2 in several fields up to 100 mT for 
H || [010] and q || [100]. With the application of magnetic 
fields, k/T continuously increases and saturates above 
H c2 . From the thermal conductivity data, we determine 
H C 2 ^40mT, in a good agreement with H C 2 determined 
from the specific heat data^ The rise of conductivity 
with field can only be attributed to the electronic contri- 
bution. As magnetic field is applied, the phonon thermal 
conductivity is either suppressed due to additional scat- 
tering from the vortices in the mixed state, or it is ap- 
proximately constant when the phonon mean free path is 
smaller than the distance between the pinned vortices. 
In contrast to the case in BaNi2As2, where the low- 
temperature k/T can be expressed as a + bT 2 in the in- 
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FIG. 3: (Color online) (a) k/T vs T of SrNi 2 P 2 in several 
fields up to 100 mT for q \\ [100] and H \\ [010]. Dashed vertical 
lines correspond to T = 0.10K and 0.25 K, at which conduc- 
tivity at different fields were determined for panel (b). (b) 
Thermal conductivity at 0.10 K and 0.25 K, divided by the 
normal state value, (reo/T)/(re n /T), of SrNi 2 P 2 as a function 
of H /H C 2- The data for fully gapped superconductors in the 
clean— (Nb) and the dirty— limit (InBi and BaNi 2 As 2 ), and 
a nodal superconductor Tl 2 Ba 2 CuOe+i (Tl-2201)— are dis- 
played for comparison. Dotted lines are guides to the eye. 



vestigatcd field region, the low-temperature k of SrNi2P2 
cannot be expressed as T a with a field-independent pa- 
rameter a. The variation of a with field may be ascribed 
to the change in the ratio of electronic and phononic con- 
ductivity with field. We therefore used the variation of 
the measured (k(T)/T)/(k„/T) at a finite low temper- 
ature as a function of field to observe the effect of the 
excited quasiparticles. 

The field dependence of the scaled thermal conductiv- 
ity of SrNi 2 P 2 at 0.1K and 0.25 K (<T C ) is shown in 
Fig. Op) as (k /T)/(kJT) vs H/H c2 for H JL q, where 
the normal state value re n is obtained at 0.1 T. Data for 
Nb^ InBi^ and Tl 2 Ba2Cu0 6 +5 (Tl-2201)^ are shown 
for comparison as representative examples of a clean s- 



4 



wave superconductor, a dirty s-wavc superconductor, and 
a <i-wave superconductor, respectively. A concave field 
dependence, clearly observed for Nb in the clean limit, 
is a consequence of states initially localized within the 
vortex cores becoming delocalized as the wave-functions 
increasingly overlap between neighboring vortices as the 
vortex density increases with magnetic field. In the dirty 
limit (InBi), where the magnitude of the superconducting 
gap is reduced by impurities, the field evolution becomes 
more gradual below i? c2 , but remains concave in the low 
field regime. In contrast, (k /T)/ (n n /T) of a nodal su- 
perconductor such as Tl-2201— has a substantial residual 
value in zero field and displays a convex field dependence 
due to the Volovik effect^ In our recent work,— we con- 
cluded that BaNi2As2 (also shown for comparison) is a 
fully gapped superconductor from the relatively fast but 
concave field dependence of kq/T in the low- field region. 
We also established that our BaNi2As2 samples were in 
the dirty limit as the estimated mean free path, Z c = 70 A, 
was much smaller than the coherence length, £ = 550 A££ 
The data for SrNi2P2 is remarkably similar to BaNi2As2. 
The value of (k /T)/(k u /T) at 0.1K is negligibly small 
(4 x 10 -3 ) in zero field, and the field dependence is con- 
cave at low fields. Thus, we conclude that SrNi 2 P2 is 
a fully gapped superconductor in the dirty limit. As in 
BaNi 2 As 2 , k /T(H) of SrNi 2 P 2 exhibits a slight shoulder 
close to H/H C 2 = 1. This could be due to a distribution of 
H C 2 in the crystal, or multiband superconductivity which 
gives rise to a similar field dependence! 39 ' 40 However, in- 
tcrband scattering in these dirty samples should wipe out 
the effects of multiband superconductivity. 

Fully gapped superconductivity in SrNi2P2 is also sug- 
gested by the comparison of the specific heat data with 
that of BaNi2As2 in zero field. Figure @] shows the 
normalized C/T with respect to normalized T/T c for 
SrNi2P2 and BaNi2As2- Both sets of data, which over- 
lap almost perfectly, can be fit with a BCS curve with 
a slightly smaller energy gap (A = 1.61 fceTc) than that 
for weak coupling (A = 1.76 /«bT c ). This is in contrast 
to the nodal superconductor Sr2Ru04^ which has a 
smaller jump at T c and, more importantly, much higher 
density of low energy excitations at low temperature. 
The data from an Fe-pnictide BaFei.84Coo.i6As2 with 
T c = 20 K— is qualitatively similar to that of the Ni- 
pnictides in Fig. 2) However, a larger specific heat jump 
indicates that a larger gap over some portions of the 
Fermi surface is required, while the larger specific heat 
at the lowest temperatures indicates that a smaller gap 
on other portions is required. Consequently, a supercon- 
ducting gap which varies significantly around the Fermi 
surface is required to explain the heat capacity data of 
Ba(Fe,Co) 2 As 2 ^— Interestingly, although it is a signif- 
icantly smaller effect, similar systematic deviations can 
be seen between the Ni-pnictide data and the single gap 
BCS fit. This indicates that weak multi-gap behavior 
might also be present in the Ni-pnictidcs, although not 
nearly as pronounced as in the Fe-pnictidcs. 
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FIG. 4: (Color online) Scaled specific heat versus tempera- 
ture for SrNi2P2, BaNi2As2,— Sr2Ru04 - an unconventional 
superconductor with comparable T c ,— and an Fe-pnictide 
BaFe1.g4Coo.i6 As2 with T C = 20K.— Solid curve represents 
a theoretical calculation based on weak coupling BCS su- 
perconductivity with a slightly smaller energy gap. The 
BaFe1.g4Coo.i6 As2 data has had the phonon contribution and 
a residual linear term subtracted. 



IV. DISCUSSIONS 

The thermal conductivity data presented here for Ni- 
pnictides display qualitatively similar behavior to that 
found in several Fc-based superconductors. In Fc- 
pnictides, a more rapid concave field dependence of the 
residual linear term than observed here has been at- 
tributed to a highly anisotropic gap.— ~— Yet other stud- 
ies of different compositions observe a finite residual lin- 
ear term with a convex field dependence, which is at- 
tributed to nodal superconductivity ! 11 ' 66 ' 67 Our data in- 
dicate that SrNi2P2 (and BaNi2As2) are fully gapped 
single (perhaps very weakly anisotropic) gap supercon- 
ductors. NQR results on LaNiAs(0,F) and Fe-pnictidcs 
are consistent with the differences observed in thermal 
transport. A Hebel-Schlichter peak is observed in 1/T\ 
in LaNiAs(0,F) followed by an exponential decay, con- 
sistent with a fully gapped s-wave superconductor 
while 1/Ti in the Fe-based materials shows evidence 
for a strongly fc-dependent ga,p£& The differences be- 
tween the Fe- and Ni-based superconductors could orig- 
inate from their contrasting Fermi surfaces, quasi- two- 
dimensional in the former, while more complicated and 
three-dimensional in the latteri 28 ' 29 ' 68 " — A recent X-ray 
spectroscopy plus DMFT study suggests that LaNiAsO 
is more strongly correlated than LaFeAsO,— although 
this seems to contradict the results of thermopower— 
and optical conductivity! 7 - 3 - From thermal transport stud- 
ies on BaNi2As2 and SrNi2P2, the superconductivity 
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in Ni-pnictidcs is most likely conventional phonon me- 
diated s-wave pairing, while the issues of the pairing 
mechanism and the gap symmetry are still controver- 
sial in Fe-pnictides. The fact that the two iso-structural 
Ni-pnictide superconductors discussed here show fully 
gapped Fermi surfaces is most noteworthy, since two out 
of their three constituents arc different. In particular, 
As and P were shown to lead to a different character of 
magnetic fluctuations in Fe-pnictide compounds^. Our 
results, however, indicate that substituting Fe with Ni 
is by far the most important aspect of the Ni-pnictide 
superconductors, leading to qualitatively different char- 
acter (and likely the origin) of superconductivity with 
fully gapped Fermi Surface. 

Finally, we want to draw attention to the high sensi- 
tivity of thermal conductivity data in SrNi2P2 to a cool- 
ing/field variation route used to reach the initial state of 
the experimental run. Figure [5] shows k(T) of SrNi2P2 
obtained for zero-field cool (ZFC) and field cool (FC) 
sample at (a) 10 mT and (b) 15 mT. For both fields the 
ZFC and FC data are markedly different. The 10 mT 
ZFC data almost coincides with the zero field data, while 
the 15 mT ZFC data is very close to the 11 mT FC 
data. Above 20 mT, the difference between the ZFC and 
FC k(T) data (not shown) is no longer detectable. Note 
that, as indicated by arrows in the inset of Fig. [5fb), the 
ZFC 15 mT k(T) deviates from FC 11 mT data around a 
characteristic temperature T-j* = 0.4 K and merges the FC 
15 mT data around T 2 * = 0.7 K (< T c ). We interpret this 
behavior in the following way: below Tf* the magnetic 
flux in the 15 mT ZFC run is frozen, and its effect on 
the heat transport is identical to that in the 11 mT FC 
sample. Above Tj" additional flux begins to penetrate 
into the sample, and the flux distribution becomes in- 
distinguishable (via thermal conductivity measurement) 
from that in the FC sample at the same field (15 mT) 
above T£. Similar differences were found for fields up 
to 20 mT although the difference between FC and ZFC 
data becomes smaller with increasing field (not shown). 
Such field hysteresis is often observed in superconduc- 
tors with strong magnetic flux pinning. Based on the 
Bean model^ which phenomcnologically describes irre- 
versible superconducting properties, in the case of zero 
field cooling magnetic flux begins to enter the sample 
from the edges, the penetration distance being propor- 
tional to the applied field in a low-field region. In the nor- 
mal state magnetic field uniformly penetrates the sam- 
ple, and therefore in a FC sample the flux distribution 
is uniform, and its thermal conductivity always repre- 
sents an upper limit for that in the ZFC sample at the 
same field. This model therefore can explain the mag- 
netic field history dependence of thermal conductivity of 
SrNi2P2 described above. However, it is also expected 
from that model that when a sample is cooled in a mag- 
netic field greater than H C 2 , and the field is then ramped 
down to the experimental field H cxp within the super- 
conducting phase, a substantial extra flux should exist in 
the sample compared to the FC case for the same field 
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FIG. 5: (Color online) k vs T of SrNi2P2 for different system 
preparation, zero-field cooling (ZFC) and field-cooling (FC), 
at (a) 10 mT and (b) 15 mT. Insets: (a) comparison of FC 
and FC* data for 10 mT, and the data at 11 mT (FC). FC* at 

10 mT data was collected after cooling the sample in a field 
above H C 2, and then ramping it down to the measurement 
field, (b) Zoomed in k vs T. 15 mT ZFC data deviates from 

11 mT data at (T*), and merges with 15 mT FC data (T 2 *), 
as indicated by the arrows. 



-f/cxp. However, the effect of such trapped flux although 
present is smaller than expected in SrNi2P2: as seen in 
the inset of Fig. 03a), the FC* data, obtained by cool- 
ing the sample in a field of 100 mT (3>iJ C 2) and then 
reducing the field to 10 mT at the lowest temperature, is 
rather close to the 10 mT FC data, and just below the 
1 1 mT FC data, indicating a difference in the energy bar- 
riers for flux to enter and leave the sample. The observed 
strong field history dependence of thermal conductivity 
emphasizes the need for field cooled experiments, espe- 
cially when one attempts to determine the symmetry of 
the superconducting gap from the field dependence of the 
low temperature properties. 
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V. CONCLUSION 

To conclude, we have performed magneto-thermal con- 
ductivity experiments on SrNi 2 P2 to identify the struc- 
ture of the superconducting gap. Both from the tem- 
perature and field dependence of thermal conductivity, 
we conclude that SrNi 2 P2 is a fully gapped supercon- 
ductor, as is the case in BaNi2As2. We stipulate that the 
fully gapped spectrum might be a universal feature in Ni- 
pnictides. In addition, we found a substantial field his- 
tory dependence of the thermal conductivity of SrNi2P2, 



which raises a note of caution for experiments done in 
magnetic field in this family of compounds. 
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